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', Stanford University 

pL^ ■ In this paper we consider probabilistic analogues of some clas- 

' sical integral geometric formulae: Weyl-Steiner tube formulae and 

, the Chern-Federer kinematic fundamental formula. The probabilis- 

tic building blocks are smooth, real-valued random fields built up 
from i.i.d. copies of centered, unit-variance smooth Gaussian fields 
on a manifold M. Specifically, we consider random fields of the form 
CLh : fp = F{yi{p),...,yk{p)) for FeC^iR^-R) and iyi,...,yk) a vector 

of i.i.d. centered, unit-variance Gaussian fields. 

The analogue of the Weyl-Steiner formula for such Gaussian- 
related fields involves a power series expansion for the Gaussian, 
rather than Lebesgue, volume of tubes: that is, power series expan- 
sions related to the marginal distribution of the field /. The formal 
expansions of the Gaussian volume of a tube are of independent ge- 
^ ■ ometric interest. 

' As in the classical Weyl-Steiner formulae, the coefficients in these 

, expansions show up in a kinematic formula for the expected Eu- 

1/^ ' ler characteristic, Xy of the excursion sets M n f^^[u, +00) — M n 

CN ; y-^{F-^[u,+oo)) of the field /. 

■ The motivation for studying the expected Euler characteristic 

comes from the well-known approximation P[suppgj^j^ /(p) > u] ~ 

5' E[x(rM^,+~))]. 

^ ! 1. Introduction. In this paper we consider perhaps the simplest non- 

S ' Gaussian models of smooth random fields on a manifold M. Our fields, 

* * I 

■ which we refer to as Gaussian related, are smooth, real-valued random fields 
built up from i.i.d. copies of centered, unit-variance smooth Gaussian fields 

' on a manifold M. Given a C^, centered, unit-variance Gaussian process y 
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on M and F G C2(M'';]R) we consider the field 

fp = F{yi{p),---,yk{p)), 

where the fields yi are i.i.d. copies of y. 
For a concrete example: set 

F{x) = \\xf, 

then the field 

f{p)=F{yi{p),...,yk{p)) 

has a Xfc marginal distribution, and has previously been referred to as a 
"x|" field [1, 20]. Note that, unlike a Gaussian field, a field with marginal 
distributions is not determined solely by its covariance function. In this work 
the fields are characterized by their covariance functions and the function 
F. 

The motivation for studying the expected Euler characteristic of the ex- 
cursions f~^[u, +oo) C M comes from the approximation 



sup f{p) > u 



E[x(r'K+oo))] 



[1, 2, 20, 22]. This approximation has found uses in medical imaging, astro- 
physics and multivariate analysis [15, 16, 20, 23]. A heuristic justification of 
the above approximation can be found in [2, 6], with a rigorous justification 
in [17]. 

Our main result. Theorem 4.1, expresses the expected Euler characteristic, 
X, of the excursion sets f~^[u, +oo) 

in terms of geometric quantities related to the curvature of the Riemannian 
structure induced by the covariance function of y [18] and geometric quan- 
tities related to the curvature of F~^[u, +oo) C M.^, viewed as a subset of 
the probability space (M'^,7jgfe), where j^k is the standard Gaussian measure 
on R'^: 

7m. (^)= / (27r)-'=/2e-l!-llV2^a;. 



A 

That is, the geometric quantities related to F~^[u, +oo) depend on both the 
curvature of F~^[u, +oo) as well as the measure 'y^k. 

In Theorem 4.1, curvature enters as coefficients in certain volume of 
tubes expansions. As these expansions are slightly nonstandard, we recall 
some basic facts about such expansions, dating back to Weyl and Steiner 
[12, 14, 19]. These expansions give an expression for the volume of a tubular 
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neighborhood of a set M C M", assumed either to be convex or an embed- 
ded submanifold, in terms of certain intrinsic measures on M, the so-called 
Lipschitz-Killing curvature (signed) measures (£j(M, •))o<j<n- These mea- 
sures depend on M and are finitely additive in M . When M is an embedded 
manifold in M", the tube formula is due to Weyl [19], and when M is a 
compact, convex domain, the expansion is generally attributed to Steiner 
(cf. [14]). The two formulae were generalized to sets of positive reach by 
Federer [11] who defined the Lipschitz-Killing curvature measures of such a 
set. 

Weyl and Steiner's formulae state that for r small enough 

Wfc(r(M, r)) ^ Hkiiy G : d{y, M) < r}) = {{y G : d{y, M) < r}) 
(1.1) 

k k 

j=0 j=0 

where d{-,M) is the standard distance function on M"; ujk = 7r'^/^/r(A;/2 + 1), 
the volume of the unit ball in M*"'; Tlk is /c-dimensional Hausdorff measure 
and Ajgfe is Lebesgue measure. It is sometimes useful to rewrite (1.1) in terms 
of Minkowski Junctionals defined by 

^-£,(M)a;,_, 

so that (1.1) reads as a (finite) Taylor series expansion 

k j 

X^,{T{M,r)) = J2Mj{Mf-. 

,7=0 

One of the deep facts about the Lipschitz-Killing curvatures, first proven 
by Weyl, is that they are intrinsic to M. That is, if we embed M into a 
different Euclidean space in an isometric fashion, the Lipschitz-Killing cur- 
vature measures of S are unchanged. They are also intrinsic in the Rieman- 
nian sense, because they are local and can be computed from a Riemannian 
metric on M. 

However, because Cj{M) are intrinsic to M, the Minkowski functionals 
of M are not intrinsic. In particular they depend on the dimension of the 
Euclidean space in which M is embedded. Strictly speaking, it is therefore 

necessary to write A^^"* (M) to clarify which Euclidean space we are talking 
about, where Xj^k refers both to which dimension M is considered to be 
embedded in as well as which measure we are using to compute the volume 
of the tube. In Section 3 we will compute the volume of a tube with measures 
other than Lebesgue; specifically we will compute a Taylor series expansion 
of the standard Gaussian volume 'j^k of certain tubes. These expansions will 
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play a key role in our analogue of the Chern-Federer kinematic fundamental 
formula (KFF) [7, 10, 11, 13] in which the Lipschitz-Killing curvatures also 
play a prominent role. 

The KFF relates the "averaged" j th Lipschitz-Killing curvature Cj {Mi n 
5-M2) to the Lipschitz-Killing curvatures of Mi and M2, averaged over "typ- 
ical" rigid motions g. Specifically, for Mi and M2, two embedded submani- 
folds of M'^, the following relation holds: 



(1.2) 



where 



Cj{MingM2)dfin{9) 



k-j 

E 

i=0 



'i+j' 




'k' 


. j . 




i 



Cj+i{Ml)£k-^{M2 



^k 



LOiLOk- 



and Gk = M'^ x 0{k) is the group of rigid motions on M*^ with Haar measure 
fik = A]Kfc X /Ifc, the product of Lebesgue measure and the invariant proba- 
bility measure on 0{k). 

The Chern-Gauss-Bonnet theorem [9, 11] states that Co{-) = x(")) the 
Euler characteristic, so that we can rewrite the case j = in (1.2) as an 
"averaged," or "expected," Euler characteristic 



(1.3) 



x{MingM2)dfikig) 



k r 

E 

i=0 



-1 



Ci{Ml)Ck-^iM2) 



keeping in mind that the measure /ifc is not a finite measure. 

This brings us to the subject studied in our work here. In this work we 
study the expected Euler characteristic of certain random sets derived from 
smooth zero-mean, unit-variance Gaussian random fields on a manifold 
M. Our main results are Gaussian analogues of (1.1) and (1.3), given in 
Corollary 3.4 and Theorem 4.1. 

Specifically, let y = (yi, . . . , yk) : M ^ M'^ be a sequence of i.i.d. zero-mean 
unit-variance Gaussian random fields (satisfying certain regularity condi- 
tions specified in Theorem 4.1 below) on a manifold M. We derive, 
using techniques similar to [1, 3, 18, 21], an expression for E[x{Mr\y-^D)], 
where D is a suitable domain in M.^. In particular, we show in Lemmas 
2.5 and 2.6 that if M is a n-manifold, with or without boundary, and 
for suitable F G C^mfc 



(1.4) 



nx{Mny-\F~'[u,+^)))]=J2Cj{M)p,iF,u), 
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where the functions Pj{F, u), for j > 1, are Gaussian analogues of A4 and 
are given by integrals of certain functions over F~^{u} = d{F~^[u, +00)) C 
M'^. In (1.4), the quantities (£j(M))o<j<n are the (total) Lipschitz-Killing 
curvature measures of the Riemannian manifold {M,g), where g is the Rie- 
mannian metric induced by the random field / 

g{Xp,Yp) = E[XpfYpf] 

for tangent vectors Xp and Yp in TpM. 

When the y^'s are defined on M" and are isotropic, the functionals Pj{F, u) 
agree with the EC densities pfj{u) oi F oy (cf. [21]), where it is shown that 
for compact domains S 

n 

(1.5) nxiS n r\u, +CX)))] = a,,nM':i^{S)p^j%f,j{u), 

3=0 

where p2 is a spectral parameter of /, namely the variance of its first par- 
tial derivatives and constants, independent of or /. It should 
be noted that (1.5) holds for any smooth isotropic field on M", not just 
Gaussian-related fields, that is, those derived from i.i.d. Gaussians. Lemmas 
2.5 and 2.6 are thus extensions of (1.5) to fields derived from i.i.d. Gaussians 
on manifolds. Further, these lemmas show that in order to compute 

E[x{Mny-\F-'[u,+oo)))], 

we need only compute (£j(M))o<j<n and pj{F,u). 

Lemmas 2.5 and 2.6 also provide a direct way to compute pj{F,u), and 
hence pfj{u) for isotropic Gaussian-related fields of the form F oy, in that 
they are represented as conditional expectations of random variables defined 
on M'^' and do not contain any information concerning the derivative of the 
fields themselves. This should be compared with earlier works (cf. [8, 20]) 
where pfj{u) (assuming without loss of generality that p2 = 1) is expressed 
as 

(1.6) /5j,,-(n)=E[det(v2/|^.J|V/|,,4 = 0,/|,,4>^/](^v/|,,(0), 

where V"^ fki,t = d'^ f{t)/dxkdxi is the Hessian of /(t), f\j is the restriction 
of / to a j-dimensional subspace of M" and ip\/ f^. ^ (0) is the density of ^ f\j^t 
for some fixed t. In these works, it was therefore necessary to work out the 
full joint distribution of the field along with its first two derivatives at a 
given point t and then carry out a fairly delicate conditioning argument. 

The connection between the functionals pj {F, u) and kinematic formulae 
is related to an extension of the Weyl-Steiner tube formulae. We discuss 
this extension in Section 3, where we derive a formal power series expression 



6 



J. E. TAYLOR 



for the standard Gaussian volume of a tube around D, that is, we derive a 
formal expression for 

(1.7) / {2.)-^/'e-\\^\\'''dX^.{x) = Y.MY{D)'-. 

The notation A^J"''' {D) is chosen to highlight the analogy between A^J**" (D), 
the coefficient of / j\ in a formal expression for the Gaussian volume of 

T{D,r) and Al/' {D), the coefficient of / j\ in a formal expression for the 
Lebesgue volume of T{D,r). Both expansions are exact in some cases, for 
instance, when D or D'^ is a compact domain. 

The "take-home message" of this work is Theorem 4.1, where we link 
Sections 2 and 3. We prove, by direct computation, that for suitable F £ 

p,{F,u) = {27r)-^/^Mf^'{F-\u,+^)), 

that is, the EC densities pjj{u) for Gaussian-related isotropic fields are (up 
to constants) coefficients in a Gaussian volume of tubes expansion around 
F~^[u, +oo). We can thus rewrite (1.4) as 

n 

(1.8) EixiM ny-'D)]=Y,C,iM)i27T)~^/^Mf'' (D). 

j=0 

We conclude with some examples in Section 5, specifically rederiving, in 
light of (1.8), the EC densities of Gaussian and fields. While these are 
not new, their derivation sheds some light on the origin of the formulae. For 
instance, although it has been shown that the EC densities for a real-valued 
Gaussian field are essentially just Hermite polynomials times the standard 
Gaussian density, (1.8) shows that the reason the Hermite polynomials ap- 
pear is the fact that they are derivatives of the Gaussian densities, hence 
(up to constants) the coefficients of powers of r in a power series expansion 
of the Gaussian measure of the tube [u — r, +oo) of radius r around [u, +oo). 

The new results in Section 5 are the EC densities of the noncentral 
fields, as well as the EC densities of what we refer to as the "correlated 
conjunction" random field which is given by taking the minimum of two 
correlated Gaussian fields defined on a manifold M. Specifically given p G 
(—1,1) and y = (?/i,y2) two independent centered, unit-variance Gaussian 
fields on M , define the new random fields 

Zl = yi, 

Z2 = p-yi + ■y2, 

Zl A Z2{p) = mm{zi{p),Z2{p)). 
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In Section 5.4 we derive the EC densities of zi A Z2 by relating the EC den- 
sities to the Gaussian measure of tubes around arbitrary cones in R^. The 
"correlated conjunctions" is a simple model for correlated random sets: the 
independent model (p = 0) has been used in certain neuroimaging applica- 
tions [24]. 

The outline of the paper is as follows: Section 2 is dedicated to deriving 
an explicit expression for the EC densities Pf,j{u) = pj{F,u). Section 3 is 
dedicated to extending the Weyl-Steiner tube formulae to measures with 
smooth densities with respect to Lebesgue measure and is possibly of in- 
dependent interest. In Section 4 we show that the functionals pj{F,u) are, 
up to constants, coefficients in a Taylor series expansion for the standard 
Gaussian measure of a tube around F^^[u, +00). Section 5 is devoted to 
examples. 

2. Euler characteristic densities. In this section we derive an explicit 
integral representation for the EC densities pj {F, u) . The proof involves some 
preliminary lemmas necessary to carry out the conditioning in (1.6). 

2.1. Gaussian random fields on manifolds. In this section we derive cer- 
tain expectations for Gaussian random fields on a manifold M. We first 
review some linear algebraic preliminaries; readers are referred to [3, 18] 
for further details. For an n-dimensional vector space V, let (T*(y),(8)) 
denote the algebra of covariant tensors on V; {/\*(y),A), the Grassman al- 



gebra and {A*{V)^ A*{V),-), the algebra of double forms A*{V)^A*{V) = 
0^,=o^'^''(^) = ®n^=o^''(^) ® ^'(^) endowed with the "double- wedge" 



The subalgebra i^'jLo^^'^ is denoted by (A*'*(^), •)• Any inner prod- 
uct ( •, •) on y induces an inner product on /\* {V) , and the trace, Tr : /\* (V) Cg" 
A*(l^)^Mis defined by 



and extended linearly. 

We call W : P) A^'^{V) a Gaussian double form if, for any basis 

By = {vi, . . . ,Vn} of V, the matrix W{vi,Vj) has (jointly) Gaussian entries. 

We recall the following useful lemma from [18]. 

Lemma 2.1. Suppose that W is a Gaussian double form; then 



product 



(a /3) • (7 5) = (a A 7) (/3 A 6). 



Tr(a0/3) = (a,/3) 



Lfc/2J 

E[iy'=] = J2 

j=0 



{k-2jy.jl2:> 



kl 
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where 

H = E[W] e A^'^V), C = E[{W - E[W]f] e A^'^iV). 

The next result we will need in subsequent sections concerns the con- 
ditional expectation of certain random double forms on finite-dimensional 
Hilbert spaces. Specifically, let {V, {■,■)) be a finite-dimensional Hilbert space 
and L CV, a subspace of V with Pl denoting orthogonal projection onto L. 
Suppose now that (^i)i<i<Af<dim(y) are i.i.d. l^-valued random variables 
with common distribution 7\/, the canonical Gaussian measure on V 

MA) = {2.^^-^^)'^ j^e~\\^\\''^dn^^^y){v). 

We evaluate the expression 

(2.1) E[a((X,,, . . . {X,„. . .,X,^))\PlXu . . .,PlXm\, 

for a G A^^^{V) and /c-tuples ik = {h, ...,ik) and jfc = (ji,. . . , j^) of {1,. . . , A^}. 

Before evaluating (2.1), we recall the definition of the annihilation oper- 
ator on T*{y), as well that of a contraction on A*{y) ® A*iy). For v , 
the annihilation operator iy on A*(^) by defining it on A^iV)^ as follows 

{iy(3){vi, . . .,Vr-l) =l3{v,Vi, . . .,Vr-l) 

and, if /3 S A^{V), we set i^,/? = 0. Although it is referred to as an annihilation 
operator, its real effect is to fix the value of the first coordinate to be v when 
the form is evaluated on r — 1 other vectors. Since A* (V) A* (V) consists 
of two "copies" of A* (V) the annihilation operator can act on either the first 
or second copy. To distinguish between these two cases we define operators 
ijv and !][, on A*{V) A*{V) by setting 

and extending linearly. 
For a£A*{V)^A*{V) 

(2.2) rl^^rl^^a = -rl^^rl^^a, 

as is easy to check from the definition. Also 

a((X,, , . . . , J, (X,, , . . . , X, J) = n '?^^.,) « 



A GAUSSIAN KINEMATIC FORMULA 



9 



so that the evaluation of (2.1) is reahy the evaluation of the conditional 
expectation of the polynomial 



1=1 



a random element of the set of linear maps from A*(V) A*{V). 

Fix {wi, . . . , ffjijiii-v^)}, an orthonormal basis for (V, (•,•)) and L C V a 
subspace with orthonormal basis {v[, . . . We define the contraction op- 
erators C,Cl,C^ on A*(y) A*(y) as follows: 



dim(y) I 



a. 



i=l 



i=l 



Cia = {C-CL)a = CL±a, 

where L-^ is the orthogonal complement of L in For any two subspaces 
Li,L2, (2.2) implies 

CliCl2 = CL2CL1 

and for a G A^'''{V) 

where a^i is the restriction of a to L. 

With the notation established, we proceed to evaluate the conditional 
expectation (2.1). 

Lemma 2.2. Suppose a e A*(^) ® A*(^), Xi, 1 <i < N, are i.i.d. 
V -valued random variables with distribution 7^, and Li, . . . , L^- are sub- 
spaces of V . Suppose further that ip = {ii, . . . , ip) and jq = (ji, - ■ ■ ,jq) are 
two p- and q-tuples of distinct elements of {1,...,A^}, arranged such that 
the first n{ip,jq) <mm{p,q) elements of ip match those of jq and there are 
no further matches in the remaining elements of ip and jq . The following 
relation then holds: 



E 



' p 



n 



(2.3) 



^n{ip,jq) N 



n vp.^^x,., 

m=n{ip,jg)+l 



n ^K^x.u. 

-m=n{ip,jq)+l ^ 
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Remark. The restriction that the indices in ip and jg be distinct is to 
rule out the trivial case when either two of the elements of ip or jg are equal 
and the polynomial 

(n.«,.)(n.^,)-. 

Proof of Lemma 2.2. We consider first the case in which p = q = 1, 
and note that, by linearity, it suffices to prove the lemma for a £ A^'IV) (8) 
A'*(y). Working from the definition of iJXi^v'xj 

n 

= {^h^'^k){Xj^,vi)a{{vk,wi, . . . ,Wr-i), {vi,w[, . . . ,w'^_^)), 

k,l=l 

where B = {vi, . . . ,Udim(y)} is some orthonormal basis for {V, { •, •)). If ii = 
ji , we can further choose the orthonormal basis for V so that {vi , . . . , Udim(Lij ) } 
is an orthonormal basis for Lj^ and {t^dim(LiJ+i) • • • j''^dim(V)} is an orthonor- 
mal basis for . The result then follows from well-known results about the 
conditional distributions of Gaussian vectors. 
Next, consider p,q> 1, 
' P \ / Q \ 



E 



f[vxj(flv'xJ»Xu...,Pu^X,,,...,PL^^X,,,...,X^ 



E 



N 



PliXi, . . .,Plj^Xm 
X E[r/x,, ri'x,^ a\X,, . . . , Pl^^X^,, . . . , Pl^^ X,,,..., Xn] 



X, 



X E 



d=2 



d=2 
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the general case then follows by iteration. □ 

Corollary 2.3. Suppose a £ A'^'^CV), Xi,...,Xf: are i.i.d. with dis- 
tribution 7y, ik and jk are as in Lemma 2.2 and vq is a unit vector in V . 
Then 

nai{X,„...,Xi^),{X,„...,X,^))\{Xi,vo),...,{Xk,vo)] 



Proof. This is just an application of Lemma 2.2 which uses the fact 
that, if the two sets ik and jk differ by more than two indices, then the 
polynomial resulting from the application of Lemma 2.2 will be identically 
zero, since the term rjy (as well as rj'y) will appear at least twice [cf. (2.2)]. 
□ 

For our purposes, we can reformulate the above corollary as follows. 

Corollary 2.4. Suppose y = {yi, . . . , yj) -.M-^W are i.i.d. mean-zero, 
unit-variance Gaussian random fields on a smooth manifold M , u is a vector 
field on W and a G A'''''{W). Then, for any (nonrandom) set of vector fields 
(Zj)i<j<m on M , where k <m< dim(M) 

E[y*a\y, {y^Zi,u),l <i<m] 

is a random double form, whose restriction to L = spanjZi, . . . , Z^} satis- 
fies, for each p £ M, 

Hy*a\L\y, {y*Zi,u), i < « < ?n] 

in the sense that for each p € M , on the set {\\i'{yp)\\ > 0} 

my*a\L{{Zh,- ■•, -^ife), (%,••• , ZjJ)p\yp, {y^Zi, u)p, l<i<m] 
- Tr-"(j'-)(a(yp)|,x(,^))/|l((Z,„ . . . , Z,J, {Zj„. . . , Z,J)^\ 

i ^ {{y*Zi,y)pf , \ 

for some universal constant Kf^j and all k-tuples {ii, ■ ■ ■ ,ik) end {ji, . . . ,jk)- 
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Proof. Without loss of generality, we can assume that the (^j)i<i<m 
are orthonormal with respect to the metric induced by any one of the y^'s. 
The result then follows from the fact that, for p fixed, (y*.^i(p))i<i<n are, 
conditional on y{p), i.i.d. Gaussian random vectors in Ty^p-^W with distribu- 
tion 'JXy^p^Rj ; combined with the previous corollary, along with the fact that 
for any vector field W on M-^, 

\'-'u^{y)'^Pv(y)Wiy)Vp^(^y^W{y)'^{y)\ < ^kj ^^(^)|p MvnW'^W ^ 

for some universal constant -fC^ j . □ 

2.2. Integral representation of pj{F,u). In this section we derive expres- 
sions for the expected Euler characteristic of the excursions of Gaussian- 
related fields, which we defined to be random fields of the form 

f{p)=Foy{p), 

where y = (yi, . . . , y^) : M — > M'^ are i.i.d. zero-mean unit- variance Gaussian 
random fields and F G C^(M'^). Recah that the EC densities were defined [21] 
for isotropic fields / on M" as follows: 

n 

(2.4) E[xiMnf-\u,+oc))]=Y,aj,nMt''",iM)f,i^'pf,,iu) 

j=0 

where the spectral parameter /_f2 is just the variance of the first-order par- 
tial derivatives of / and pfj{u) depend only on the finite-dimensional dis- 
tributions of / and its first two derivatives at any spatial location. Suppose 
fi = Foyi and f2 = Foy2 are two (not identically distributed) isotropic 
Gaussian-related fields with 1x2,1 = /^2,2, that is, with equal variance of the 
first derivatives. Then, by inspecting the results in [20], it can be shown that 
the expected Euler characteristics, /~^[n, -|-oo), will be identical (this will 
in fact be proven in the next lemma). In particular, the EC densities are 
functionals of F, so that, for isotropic Gaussian-related fields we can rewrite 

(2.4) as 

n 

(2.5) ¥.[x{M n r\u,+^))] = Y,a.j,nMn-j{M)iJj^p,{F,u) 

for some functionals pj, < j < 00. 

The simple case k = 1 and F{x) = x or F = ld, the identity map, leads to 
a real-valued Gaussian field, for which it is known [1] that 

Pj{ld,u) = (2^)-(j+i)/2 r Hj{x)e-'''/^dx, 

J u 
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with Hj the jth Hermite polynomial. In the manifold setting, if / is a real- 
valued zero- mean unit- variance Gaussian field on a manifold M, it was 
shown in [18] that 

n 

E[x(M n f-^[u, +00))] = ^ q{M)pjiid, u), 

j=0 

that is, that (2.5) holds with aj^nMn-j{M)V /"^ replaced by Cj{M), where 
the Riemannian metric with respect to which the Cj{M) are computed is 
the one induced by /, that is, 

which is assumed to be C^. 

One way to think of the relation (2.5) is as a partitioning of the geometric 
information from the parameter space and information depending on the 
distribution of the field. The result in the real-valued Gaussian case suggests 
that this partitioning might hold for Gaussian-related fields on manifolds. 
This would mean that in order to compute expected Euler characteristics, 
we would simply need to calculate the Lipschitz-Killing curvatures of (M, g) 
(recall that g is the metric induced by any one of the coordinate random 
fields yi) and the EC densities of a corresponding Gaussian-related isotropic 
field. The following lemma shows that this is indeed the case, and our work is 
reduced to computing the EC densities of isotropic Gaussian-related fields. 

We set VfE{p) = (Eifip), . . . , Enfip)) for some section of ©(M), the 
bundle of orthonormal frames on (M, g) and 

ae{y) = (w„e")"^iB^„(o,e)(y) A '^yj^ ' 

an approximation to the Dirac delta, as in [18]. 

Lemma 2.5. Given y = (yi, . . . ,yk) i.i.d. suitably regular centered, unit- 
variance Gaussian random fields on M, a n-manifold with or without 
boundary, and F S C'^{M.^) is such that f = F o y satisfies: 

(i) lim,^oIE[/A/V/|;(a,)l|j>„|] = /^,lim,^oIE[V/|;(a,)l|j>„}], 

(ii) every critical point of f in f^^{u — 6,u + 6) for some 6 > is non- 
degenerate (cf. [18] J, P-a.s. 

Then 

n 

E[x(A/n rHn,+oo))] =Y,Cj{M)pj{F,u) 

j=0 

for some functionals pj (F, u) . 
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Proof. We will prove the statement when M has no boundary, omit- 
ting the calculations for M with boundary, as they are similar to those in 
Theorem 5.1 of [18]. Further details can be found in [3]. 

Following the arguments in Theorem 4.1 of [18], under the assumptions 
of the theorem, 

E[x(MnrH^,+^))] 

limE[V/^(a,)l|^>„}] 
= /,,S^[l{/>".llv/.il<^}Tr"'((-V'/)")] VoIm,. 

= - f limE[l|^>„,|,v/,|,<,}E[Tr^^((-V2/)")|/,V/E]]VolM. 

We will therefore devote our attention to evaluating 
E[(-V2/n/,V/ij]EA"''^(M), 
which we will show, for each p £ M is of the form 

(2.6) E[{-V'fr\f,VfE]p= E -:RiMp) 



9 ■ 



j=0 



for some random double forms Aj{p) [measurable with respect to cr{f{p), 
V/£;(p))]. Since the distribution of the random vector {f {p) f e{p)) G M""^-*^ 
is independent of spatial location, this will prove the claim. To see that the 
distribution of the pair (/(p), V/£;(p)) is independent of spatial location, 
note that 

f{p)=F{yi{p),...,yk{p)) 

and 



^ dF 



yfE,iip)=E^f{p)=Y.E^yJ{p) . 

i=i 



(yi(p)v,?/fc(p)) 



Since E is an orthonormal set of frames, the Eiyj[p) are i.i.d. standard 
normal random variables, as are the iyi{p), ■ ■ ■ ,yk{p)) and the two sets of 
Gaussian random variables are independent of each other. 

We now turn to (2.6). Conditional on y and y*, V^/ is a Gaussian double 
form. Further (cf. Section 2.4 of [18]), 

E[V'f\y,y*]=y*V'F-l{y2y, 



dyi J 

y*V^F + I{VF,-V\\xf/2)\,-_ 
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= y-V^F + I{VF{y),y), 



E[(V2/ - E[V'f\y, yS'\y, y,] = -{I^ + 2R) ^ 

= -{I^ + 2R)\\VFf. 



By Lemma 2.1, it follows that 

i-E[(-v2/ri2/,y1 

W2J (_-\\n+k 

= E (n - 2k)im^ ^y*^'^ + {VFiy),y)ir''\\VFni' + 2R) 

(2-7) = E 7t(-^)'II^^II" 

j=o -^^ 

/=o 

^ J3^._;^, ^n-2,-^((VF(y)/||VF(^)||,y))J"-^^- 

By Corollary 2.4, 

E[(y*v2F)'|y, V/^;] = E[(y*v2F)'|y, {y,Ei,VF), . . . , (y,i?„, VF)] 

= Tr^^^(V2F|vi.x)'/' + ErrKV/i=;,y), 

where 

A'''(M)9ErrKV/s,y) = 0(||V/s|||„||(V2F)'|U2.^J 
in the sense described in Corollary 2.4. Combining this with (2.7) 
1 



n! 



]E[(-vV)"|y,v/i,] 

Ln/2J 



^ (-i2)^/-2i||v^||n 
V .7=0 



X iTV'^^-V^fivF^/IIVFII)'") +Eit"(V/e,!,) 
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where, up to constant factors 

Ln/2J n-2j 

Err"(V/£;,y)= ^ ^ ^^^^^^^^^^^^ 

j=0 1=0 

from which the claim readily follows. □ 

We have now reduced the problem of calculating expected Euler char- 
acteristics for Gaussian-related fields to the isotropic case. The following 
lemma gives an explicit integral representation for these EC densities. 

Lemma 2.6. Suppose f = Foz satisfies the conditions of Lemma 2.5 
where z = {zi, . . . , Zk) are i.i.d. isotropic, zero-mean unit-variance suitably 
regular Gaussian fields on M" that induce the standard Riemannian structure 
on M". Further, suppose that, for some 6 > 0: 

(i) there exist Ci , C2 > such that on F~^{u — 6,u -i- 5) 

Ci<\\VF{x)\\<C2; 

(ii) \/F is Lipschitz on F^^{u — 6,u + 5); 

(iii) the functions 

Gn,l,F{z) 

= I ||VF(x)||G„,^,^(x)(27^)-'=/2e-li-•|l'/2rf?^fc_,(^) 

\\VF{x)\\Hr,.i-i{{VF{y),y){x)) 
xTr^^^(^)(-V2F(x)/||VF(2;)||)'(27r)-^/2g-||x||V2^-^^_^(^) 

are continuous on [u — 5,u + 5) ; 

(iv) for all n, I 

Jim ^E[l||^(j^)_„l<,||i/„_i_K-(VF(y),y))||V2F(x)|U2,j^,|] < 00. 
Then, 

Proof. Using the original representation of the EC densities in [1] for 
random fields on M" 

Pn{F,u) = pf,,{u) 
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lim 



^ ) det(- v2/|„_i) f = u, V/|„_i = 



1 



-E 



L{|/-«|<£}l{||V/|„_i||<e}(^^ 



V'/,V/|„_i('",0) 

^det(V2/|n-i) 



where 



df 



df 



is the vector made up of the first n — 1 coordinates of V/, '/^/,v/|„_i is the 
joint density of (/,V/|„_i) and V^/|„_i is the matrix of the Hessian of 
/ in the first n — 1 coordinates. We can rewrite the expression inside the 
expectation above as 

f)'cle.(-VV,..-.) 



1 



;{{VF,y,{d/dtn))r 



{n-iy. 

X i-v^fr^Hid/dh,. . . {d/dh,. . .,d/dtn-i)). 

From calculations similar to those in Lemma 2.5 (just restricting V^/ to the 
first n — 1 coordinates and reapplying Corollary 2.4), we see that 



E 



Ifj det(-VV|„-0.,^,l<^<n-l 



df 



iivi^ir 

(2^0^ 



1=0 



I 



((VF(y),y)) 



xTV^^ (-V^F/||VF||)'+Err"-^(V/|„_i,y) 

n.-l 

,n-l/ 



n — 1 
/ 



G„,z,i.(2/)+Tr(Err"-^(V/|„_i,y)), 



where 



Err"-^(V/|„_i,y) 

n-l 

= 5]/-i-'i?„_i_K(VF(y),y))O(l|v/|„_i||^„-i||(v2F)'||«.,K0- 

Assumptions (i) and (iv) above imply that 
1 



lim 

£-+0 2e'"a;„_i 



E[l||^(,)_„|<,|l|||;i_^,|<,}|Tr(Err"-i(V/|„_i,2/))|]=0. 
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Passing to the limit for the remaining terms 
PF,n{u) = hnr 



71-1 

x^(-ir"i-' 



«=0 



n — 1 



X E 



E 



{||V/|„_i||<e}l{|/-«|<e} 



l|VF(y^ 
(27r)i/2 



■Gn,l,F{y) 



1 



n-1 



:hm — Vf-l) 



n — 1 



where 



X E 



7ffin-i(5K„-i(0,e/||VF(y)||)) 



^„_i(./||VF(y)||)- 



By assumption (i), 



lim D{y, e) = {2Tr) 

e— >0 



-(n-l)/2 



for all y G F ^(n — (5, ti + (5). Further for all such y, 



|D(y,e)-(27r)-("-i)/2|< sup 

0<r<e/Ci 



7R"-i(-BMn-i(0,r)) 



„n-l 



-(2vr)-i' 



(n-l)/2 



where the function L{e) is increasing and continuous on [0,+oo) and 



L 



limL(e) = 0. 



Therefore, for e < 5, 



-E[l||^(,)„„l<,|D(y,£)||VF(2/)||G„,;,,.(y)] 



(2vr)- 



-(n-l)/2j_ 



2e 



E[l||^(j^)_„l<,|||VF(y)||G„,,,i.(y)] 



< 



L(£: 

2e 



•E[l||^(j^)_„l<,|||VF(y)|||G„,,^(y)|]. 



The above implies that if, for some 5 > 0, 

(2.8) ^^vJrn^{\F{y)-u\<e}\^ny)\\Gn,lAy)\ < 

0<e<5 
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we need only evaluate 

limlE[l||^(,)„„l<,|||VF(y)||G„,,,i.(y)]. 

The Lipschitz assumption on VF allows us to apply Federer's coarea 
formula [5, 11] so that 

lim^E[l||^(y)_,|<,}||VF(y)||G„,,,^(y)] 

= li^iTT / Gn,l,F{z)dz 
= Gn,l,F{u). 

The Lipschitz assumption on V-F and the coarea formula also imply (2.5). 
□ 

3. Steiner formulae for Lebesgue and other measures. In the previous 
section we derived an integral representation for pj{F,u). Ultimately, we 
will connect this integral representation with a certain volume of tubes ex- 
pansion, to which this section is devoted. In this section we generalize the 
Weyl-Steiner tube formulae in order to compute non-Lebesgue volume of 
tubes, in particular Gaussian volumes. 

We will limit ourselves to tubes around domains D in M.^ with outward 
pointing unit normal vector field v. That is, D is a closed set with nonempty 
interior such that dD is an embedded hypersurface. We write 

T{D, r) = {y G : d{y, D) < r] 

and 

Pd (y) = arg min d{x, y) , 

the metric projection onto D when the minimum is unique. We recall the 
definition of the critical radius of D, 

rc{D) = sup{r > : d{y, D) <r ^ Poiv) is unique}, 

and we will assume throughout that rc{D) > 0. 

The following lemma, whose proof we omit, summarizes some facts about 
the distance function of D, particularly the growth of the eigenvalues of the 
shape operator of the hypersurfaces at distance d from D. For a hypersurface 
M cM.^ with outward pointing normal u and y ^ M the shape operator Sy 
at y is defined by 

Sy{Xy,Zy) = - {X y , V Z y 1^) , 

where V represents the usual covariant derivative on M'^. 
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Lemma 3.1. Suppose F£){y) = d{y,dD) for some domain D. Then, 

where Tx :T*{TyM.'') — > T* (Tx+y^^) represents translation from one tangent 
space to another. Further, let \i{y), . . . , Xn-i{y) be the eigenvalues of the 
shape operator Sy acting on TyF~^{d{y, D)} ; then 



l-d(i,,D)\i(PD{y)y 



Lemma 3.1 describes the geometry of the hypersurfaces at distance r 
from D and can be used to derive Steiner's formula, Theorem 3.2, which 
gives an expression for the Lebesgue volume of the tube T{D,r) around D, 
in terms of the Lipschitz-Killing curvature measures •))o<j<fc of D. 

The interested reader is referred to [11, 12, 13, 16, 19] for additional details. 



Theorem 3.2 (Steiner's formula). If we set 

Ar = {{p, Xp) G N{dD) :peAndD,Xp = rup}, 

where N{dD) is the normal bundle of dD in M'^, then Steiner's formula 
states that for r < rc{D) 



(3.1) 



Wfc_i(expa^(A)) = T—^M^ {D, A) 



In the next corollary, which is just an application of Steiner's formula, we 
derive a formal Taylor series expansion of the integral of smooth functions 
over T{D,r). The functions are assumed to be in 5(]R'^), the Schwartz space 
of functions on R*^. 



Corollary 3.3. Suppose D is a domain in M'^ with critical radius 
rc{D) > 0. Further, suppose D satisfies, for any <j <k. 



I T-7rM-^J*'(Ad^)<oo, 
JdD 1 + IfII'^ ■' 



for some /3 > 0. Then, for any n, r < min(rc(-D), l),(p G 5(]R'^) there exist 
{Ci{D;(p))i>Q and a constant Kn{D,ip) such that 



(3.2) 



I ^dX^k- f ipdX^k -^^Ci{D;ip) 

JT{D,r) JD ^ t! 



< 



'n + 2)l 
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Specifically, 



i-i 



Ci{D;^)=J2 



m=Q 



dD 



1\ 



1—m 



Mtt,{D,dx), 



where d/dv represents differentiation in the direction of the outward unit 
normal v and {D, •) is defined to be zero for j > k. 

Proof. As (/? G 5(M'^') for n > 0, at each p G dD, ip{ex.pQ^{p,rVp)) can 
be expressed as 



(3.3) v'(expgo(p, rvp)) = 



(n + 1)! dz^"+i 



o{r,p) 



where a{r,p) = p + 9{r,p)rup for some 6 : [0, r^D)] x dD (0, 1). 
Combining (3.1) and (3.3) we see 



expso(Ar-) 

n+k~l I I 



Lp{y) dUk-iiy) 



71 

z=o m=0 
k-l 



I 



A du"" 



M'-l^,{D,dp) 



1=0 



,n+l 



and 



fc-i 

i=0 



A (n+ 1)! dz^'^+i 



a(r,p) 



(n + 1)! 



o(r,p) 



„n+l 



< 



(n + 1)! i<i<k 

k-l 



< 



X 

1=0 



Ck max sup( (1 + 
1 



n+1 



l + |a(r,p)|/3 
K„,(D,(/p). 



(n + 1)!' 

Integrating over dD and [0,r] we are left with (3.2), which concludes the 
proof. □ 

The following corollary gives the expression for Ci{D;ip) when (p is the 
density of with respect to Aj^fe, that is, 

= (^' 



HI^IIV2 
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The Ci{D;ip) are coefficients in the promised power series expansion of the 
Gaussian volume of the tube T{D,r). 

Corollary 3.4. //tk* is the canonical Gaussian measure on M.^, we 
define, for 1>1, 

X H^i{p,u,))e-MV2M^^^jD^dp) 

and 

If further, 

for some (3 > 0, then the following formal expansion holds: 

OO j 

7k. {TiD, r)) = 7iK. {D) + ^ ^-M]^' (D). 

j=i J- 

3.1. Normalization o/A^J**'(-). We conclude this section by showing 
that the functionals A^J*'°(-) can be seen to be normalized independent of 
the dimension k, and can thus be extended, at least formally, to sets in M^. 
Because of this normalization, we would be justified in writing M'^ instead 
of M^v^*" in what follows, though we keep the latter notation in order to 
minimize confusion. 

For a finite subset I of N, we denote the projection from onto M"^ 
by TTj. For D C M*=, with k = #T, we can define the functionals A^J(-) on 
cylindrical sets of the form ttj^{D) as 

M]i7r^\D))^Mf'iD). 

The functionals can then hypothetically be extended by taking limits 

of cylindrical sets, though we will not pursue this matter here, except to show 
that these limits exist for at least j = 0, 1, 2. Of course, the functional Mq is 
well defined, since it is just the infinite-dimensional i.i.d. Gaussian measure 
on M^. 

For j > 0, however, it is not immediately clear that the extension is possi- 
ble. The following corollary, which is really just an application of the coarea 
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formula in [5] and whose proof we omit, shows why such extensions are pos- 
sible, at least in the cases j = 1,2. By the coarea formula the expression 
for A^2"'' (-^) is the integral of what is referred to as the mean Gaussian 
curvature of dD in [4], and can be extended to codimension-1 submanifolds 
of Wiener space and other infinite-dimensional Gaussian measure spaces. 



Corollary 3.5. Suppose the domain D is given by F ^[u,+oo). If 
both MI^\F-\u, +oo)) and7W^«'(F-i[n,+cx))) <oo, then 



X7'(F-i[u,+cx))) =E[||VF|| \ F = u]^f{u), 



MY {F-'^[u,+oo))=¥. 



-LF + 



V'^F{VF,VF) 

Wfv 



F = u 



(pFiu), 



where (fpi-) is the density of the random variable F{Zi, . . . , Z^) for Z,, i.i.d. 
N(0,1) and 

is the Ornstein-Uhlenbeck operator on C^(]R'^). 



4. EC densities and Weyl Steiner formulae: the Gaussian KFF. In Sec- 
tion 2 we derived an expression for the EC densities of Gaussian-related fields 
involving the integral of certain functions Gn,i,F over the surface F~^{u}. 
In Section 3 we derived expressions for coefficients in a formal series expan- 
sion of the canonical Gaussian volume of a tube around a domain D. 
In this section we will show that these expressions agree, up to a constant. 
Specifically, we show that for suitable F G C^(M^) 

which by Lemma 2.5 translates into the following Gaussian KFF: 

» n 

(4.1) / x(M n y{u)-'D) d¥{u) = Y,{2t:)-'J'^C,{M)M'Y (D). 



Theorem 4.1 (Gaussian KFF). Suppose y = {yi, . . . ,yk) are i.i.d. real- 
valued zero-mean unit-variance suitably regular (cf. [18] j Gaussian random 
fields on a n-dimensional manifold M that satisfy 



sup \\yi{p) 

q&Br{p,h) 



2>e 



o{h'' 
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for any h < Hq and any e > where 

Up) = {y^{p)yy^,E{p)y^yi,E{p)) g iR x m" m-, 

\\{X,V,H)\\2 = \X\ + ||y||Rn + ||iJ||552jg,„, 

and VfE (resp. V'^fE) are the coefficients ofVf (resp. V^/J read off in a 
fixed orthonormal frame {Ei ,...,£'„). 

Let fo = Fe) o y where Fd is the distance function of a domain in Mf^ , 
not necessarily compact, with rc{D) > which satisfies: 

(i) the functions 

are continuous on some neighborhood of zero where v is the outward pointing 
unit normal vector field; 

(ii) the second fundamental form S of dD is hounded, that is, 

\S{X,Y),\<K\\XMYA 
for some K > and all x e dD. 
Then, 

E[x{M n [0, +oo))] = E[x(M n y-'D)] 

In particular, the above relation holds for every compact domain D, 
and its complement D^. 

Remark. The conditions above are not necessary; indeed, there are 
cases such as the F field defined in [20] (cf. Section 5.3) where the bound- 
edness condition on the second fundamental form of the boundary of the 
domain is not satisfied but the EC densities exist and are given by 

n-i / _ 1 \ 1 

E(-ir"'^' [ I j lmi-E[l{|^(,)„,|}||VF||G„,,^(y)]. 

However, it should be noted that in this case, the above coefficients are 
not the true coefficients in the corresponding power series, because the do- 
main Dy_, defined in Section 5.3, has rc{Du) = 0. 
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Remark. Using the generalized Morse theorem of [16] or stratified Morse 
theory [3], Theorem 4.1 can be extended to include piecewise smooth do- 
mains and/or submanifolds of M'^ (i-e., the domains in the space where the 
random fields take values). We give an example of a piecewise smooth con- 
vex domain D in Section 5.4, where we calculate the AlJ'*'° (D) for D a cone 
in M? and show that the result agrees with known results for right-angled 
cones [24]. One approach to these generalizations, following the notions of 
continuity of the Lipschitz-Killing curvature measures as in [11], would be 
to construct a limiting argument to justify the following computation for 
certain smooth D sets with positive critical radius in M.^: 

Jn 

= lim / xiMr\y{uj)-'^T{D,r))dF{uj) 

In the interest of brevity, we will not pursue these generalizations here. The 
above theorem, however, does specify the functional form that all of these 
generalizations should have, that is, the contribution of the parameter space 
is in the form of the Lipschitz-Killing curvature measures, which can be 
defined for piecewise submanifolds of manifolds, and the contribution 
from the Gaussian space is in the form of the coefficients in an expansion for 
the volume of a Gaussian tube, which can similarly be defined for piecewise 
manifolds in M'^. For a more geometric approach to the above problem, 
see [3] which shows that (4.2) can be thought of as a limit of the classical 
KFF. 

Remark. Note that the conditions on the y^'s are not overly restrictive 
(cf. [18]), as fields whose second derivatives have a covariance function 
satisfying the "usual" l/(— log(/i))^^'^ moduli of continuity conditions are 
included above. 

Proof. Most of the work has been done in the preliminary lemmas. 
The conclusion 

Pj{FD,u) = {27:ril^Mf\D) 

follows from the fact that Jd satisfies the conditions of Lemmas 2.5 and 2.6. 
Lemma 3.1 ensures that the second fundamental form, or shape operator 5 
of dD is bounded. Verifying the equality above then follows from compar- 
ing the definitions in Corollary 3.4 of jvC-^^ [D) and those of Gn^i^Fo (0) 
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Lemma 2.6, noting that for x G dD and all Xx,Yx G T^dD, 

VFoix) = u^, 

The fact that the conditions of Lemma 2.5 are satisfied follows from Lem- 
mas 2.4 and 2.5 of [18] combined with the growth conditions on yi and 
its derivatives as well as the boundedness of S. As for the conditions of 
Lemma 2.6, since is a distance function of a set with positive critical 
radius, conditions (ii) and (iv) (cf. [11]) of Lemma 2.6 are automatically 
satisfied. A straightforward calculation shows that assumptions (i) and (ii) 
above imply conditions (ii) and (iii) of Lemma 2.6. □ 

5. Examples. Throughout this section, y = (yi, . . . , yk) : M — > M''" will de- 
note a generic sequence of i.i.d. zero-mean unit-variance Gaussian fields on 
M satisfying the conditions of Theorem 4.1, whose dimension k will vary as 
needed in each example. 

5.1. Real-valued Gaussian processes. Given a unit vector z £M.^, we de- 
fine the function tpz{x) = {x,z), so that o y is a real-valued, centered 
unit-variance Gaussian random field. As mentioned in Section 2, the EC 
densities of ■0^ ° y are given, for j >1, by 

(5.1) M^z,u) = ^-l^H,^,{u)e-^'/\ 

This result can easily be rederived in light of Theorem 4.1 as follows. Let 
Du = ijj^^ {[u, +00)) be a half-space in M'^. Clearly 

T{Du,r) = ip~^[u - r, +00) = D^-r, 

that is, a tube around a half-space in is another half-space. As 

^e-^'/' = i-iyH,ix)e-y\ 

it follows that 

j^,{T{D^,r)) = l-^u-r) 



00 j 1 



l-$(n) + E-^F,_i(n)e- /^ 



so that, for j > 1, 



1 



Mf\D^) = ^U,^^{u)e- 
which, by Theorem 4.1, implies (5.1). 



v?j2 
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5.2. The o.'^d noncentral case. In this subsection we derive the 
EC densities of a xt random field, as defined in [1, 20]. Note that the EC 
densities were also derived in [20], but we rederive them here as a simple 
application of Theorem 4.1. We set G{x) = then g = G o y is a, xl 

random field on M. Note that the EC densities of ^ are related to the EC 
densities of g by 



Py^jiVu) = PjiVG, y/u) = pj{G,u) = pgj{u), 

so that it suffices to calculate the EC densities of y^. 
We set 

-1, 



D^ = VG [x,+oo)=R>'\B^k{0,x), 
so that T{Dx,r) =Dx-r- Therefore, 

7ufe(T(L>^,r)) = 7Kfc(^)x-r)• 
It remains to express the right-hand side above as a Taylor series in r. 

The density fk = d{^/G^{'^^k)) / dX^. of the square root of a x\ random 
variable is 

' r(yfc/2)2(^-2)/2^ 

and 

oo 



t=x 



Direct calculations show that 

dP-^ r(A:/2)2('=-2)/2 ^ l{fc>i-m-2Z} y^_i_^_2l 



1=0 m=0 



Combining the two equations gives 

7Rfe (D^-r) = TRfe {Dx) + E -[r(fc/2)2(^-2)/2 



4(i-i)/2Ji-i-2/ 

E E ^{k>j-m-2l} 
V i=0 m=0 



j — 1 — m — 21 



m!/!2' 



X — ^ _^2m+2l 
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We thus conclude, by Theorem 4.1, that for j > 1, 



f.k jp 12 



PjWG,x) (27r)i/2r(fe/2)2('=-2)/2 

{k>j-m-2l} [^j_i_^_2l 



4{i-l)/2Ji-l-2/ 



\ 1=0 m=0 

and thus the EC densities for j > 1 of the xi I'andom field g are given by 

if(fc-j)/2e-V2 



/L(j-i)/2Ji-i-2i 



1=0 m=0 



j — 1 — m — 21 

+1 

m\lV2} 



which agrees with [20]. 

Using the above formula for the EC densities, we can derive the EC 
densities of a noncentral \^ which we define to be a Gaussian-related field 
with the function G^, : M'' ^ M 

Gfiiv) = \\y- ^lf. 

Recalling that the density fa^k of the square root of a noncentral random 
variable with noncentrality parameter a can be expressed as 

OO 

fa,k{x) = -/fe+2j(x), 

i=o ^■ 

where fk{x) is as above, the density of a xt random variable, and noting 
that, in our case a = we obtain the following. 



Lemma 5.1. 



Q 2H\ (27r)i/2r((A; + 2i)/2)2(^^+2»"2)/ 

/L(i-i)/2j j-l-2« 

^{k>j-m-2l-2i} 



1=0 m=0 



k + 2i-l 
j — 1 — m — 21 



(-iy-i+-+'(j-l)! ..^+; 

X. ~~ — j (J) 

m!n2' 
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5.3. The F case. In this subsection we derive the EC densities of an 
FkiM random field, as defined in [20]. Unhke the previous two subsections, 
instead of using the representation of the EC densities as coefficients in a 
Taylor series expansion of the volume of certain tubular neighborhoods, for 
the Ffci,fc2 fisld we use the expression for mY'^"^ (•) given in Corollary 3.4. 

We set 

k2 yUlm? 



F{y) 



^1 Ttiivk^+iY' 



and define the -Ffci^fcj random fields f hy f = Foy. Setting = F ^[u, +oo) 
we see 



dDu = F-'{u} = U Srk.u/k.O^"') X Sr 



reM+ 

The definition of the functionals (A^J**"^^*"^ {'))j>Oi along with the fact {\7F{y), 
y) = for all y £ R'' and 

{0, n odd, 

imply that 
{27ry/^p,{F,u)=M]^'''^'^Du) 

= (J-1)! E ^ ^ 



=0 

X 



(5.2) 



112'- (27r)('=i+'=2)/2 
JdD„ 



L(i-i)/2j , ..I . 



1=0 



/!2' (27r)('=i+'=2)/2 



X / e "-11^/2 



2/9 1 



IdDu U- 21 -1)1 

xTV^^"(5^^2j-i)dH,,+fc,_i(x). 
We now proceed to evaluate, for 0<m,<A:i + A;2 — 1, 

Tr^^-(5g),J(y)=Tr^^-(^)((-||VF(y)||-iv2F(y)|eD,(,,r) 

in terms of U{y) = EjLiiVjf and V{y) = EjliiVj+k,? and G{y) = U{y)/ 
V{y). 
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Lemma 5.2. We have the following expression for Tr^^^ (Sq^^^) 



m! ' m! 



m — I I \ I 



.^VG{l + G) 
Further, 

_ r{{ki + k2-m-i)/2) fcm\-(^^+^^-^)/^ 



xE(-i) 

i=0 



k2 ) \m - i ) \ i 



Proof. Since F = k2U/kiV, a straightforward calculation shows 



\VF\\-^VF=J ^ ^ ^^ VU-J ^^ ^V, 



1 



2v/G(l + G) 

X (^V^C/ - -^(dy ^dU + dU^ dV) 



2G G 
+ -^dV(^dV-^V^v]. 

y3/2 /y 



We now evaluate the matrix of ||VF||~^V^-F in a specific set of orthonor- 
mal frames of VF-*- . Considering R'^i+^^z i-,g the product M'^i x M'^^ ^[^i^ the 
product metric, we fix subspaces Li the kernel of VU in R'"'^ and similarly L2, 
the kernel of in M.^^, for which we choose orthonormal frames Bi,B2. 
The desired set of frames for the kernel of VF is then B = {Bi,B2,X}, 
where 



The matrix of — 1| V-F||^^V^F in this set of frames is diagonal, with entries 

k-^ — Itimes ^2 ~ Itimes 



1 



i-l,...,-l,G,...,G,o), 
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from which the result follows by expanding the trace of the mth power of 
such a diagonal matrix. 

Federer's coarea formula (cf. [5, 11]) implies 



E 



|VF|| — T¥^^^(5: 



dDf 



F = u 



-{u). 



The second conclusion now follows after substituting in the density 
dF*{lM.H+k2) / dXwL and noting that V{1 + G) = U + V ^ x\-^+k2 '^^ indepen- 
dent of G and 



r((fci + fc2)/2+p) 

r((A;i + A;2)/2) ' 



□ 



Combining (5.2) and the previous lemma we have, in agreement with [20], 
the following 

Lemma 5.3. The EC densities for the Fk^^k2 random field are given by 
T{{ki + k2-j)/2) /fciuV^i-^')/^/ kiuY'^^^+^^-'^y'^ 



20--2)/2r(A;i/2)r(A;2/2) V k2 



1 + 



^ ' Ti{k, + k,-j)/2)ll 



i=Q 



ki - 1 

j - 1 - 2/ - i 



A:2- 1 



5.4. Cones in M? and correlated conjunctions. In this section we study 
conjunctions of correlated Gaussian fields, which in terms of fields can be de- 
fined in terms of the minimum of two correlated Gaussian processes. Specif- 
ically, as in the Introduction, given y = (1/1,2/2) i-i-d. centered unit-variance 
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Gaussian processes on some manifold M, we form two new Gaussian pro- 
cesses as follows: 

zi = yi, 



and define zi and Z2, our isotropic versions of these processes, as in Section 3. 
We define the conjunction of zi and Z2 at a point p G M by 

ziAz2ip)=mm{zi{p),Z2ip)). 

It is easy to see that 

zi A Z2^[u, +oo) = Zi^ [u, +oo) n Z2^[u, +oo) 
= y~\K{u,p)), 

where K{u,p) is a cone in so that, to calculate the EC densities for 
this process, it suffices to calculate the expected Euler characteristic of 
M n y^^K for a general cone with arbitrary vertex in M?, which we pro- 
ceed to do. 

Define the cone 

C{vi,V2,w) = {z : z = w + aivi + 02^2, «i, ^2 > 0} 

= w + {z e TyjR^ : z = aiVi + 02^2,01,02 > 0}, 

where T^M^ is the tangent space to at w. A sketch of the cone appears 
in Figure 1. We derive an expression for the quantity 

M^Civi,V2,w)) 

which, by Theorem 4.1 and the remarks following it, allows us to compute 

K[xiMny-^C{vuV2,w))]. 
Associated to C{vi,V2,w) is its normal cone 

C-^{vi,V2,w) = {zeT^R^:{z,vi) <0,{z,V2) <0}, 

with link 

L{C^{vi,V2,w)) = {ze S{T.^M.^) : {z, vi) < 0, (z, V2) < 0}, 

where S'(r^R^) is the unit sphere in T^M^. 

In this case, the domain C{vi,V2,w) is not smooth; however, for any 5 > 0, 
T{C{vi,V2,w),S) is at least and the coefficients of 

7r2 {T{T{C{vi,V2,w),d),r)), 
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• (0, 0) 

Fig. 1. Tube of radius r around C{vi,V2,w) . 



as a power series in r will, for small 5 > 0, be close to those oi^^2(T{C{vi,V2,w), 
r)), that is, for small 6 > 0, 

E[x(M n y-'C{vi,V2,w))] ~ E[x{M n y-\T{C{v,,V2, w),S)))] 

1 " 

\ ^) j=o 



where, for the cone C{vi,V2,w), M^^^ {C{vi,V2,w)) is defined as the coeffi- 
cient of r-'/j! in a Taylor series expansion of 'y^2{T{D,r)). The expansion is 
split up into integrals over three regions, depicted in Figure 2. 

Without loss of generality, we can choose a basis of T^jM^ so that the coeffi- 
cients of vi are (1, 0) and those of V2 are {cosO, sinO) where 9 = cos~^((?;i, ^2)). 
We then set Vi to be the unit vector orthogonal to vi such that {vi,V2) > 0, 



34 



J. E. TAYLOR 




• (0, 0) 

Fig. 2. Regions of integration for power series expansion of j^2{T(C{vi,V2,'w))). 



that is, the coefficients of Vi with respect to our chosen basis are (0, 1). The 
volume of the first region, Ri{r), is thus 

7K2(i?i(r)) =7iR2({z G R'^:{vi,z) G [{vi,w) - r, {v^ ,w)], {vi, z) > {vi,w)}) 



2vr V.tl J 



By symmetry, 

lR2{R3{r)) =-f^2{{z e R'^:{v2,z) G [{v2,w) - r, (u^, w)], (?;2,z) > {v2,w)}) 



2vr V.tlj! 
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Finally, using polar coordinates 



2vr ^J[0,r]xL(Ci{t;i,i)2,»«)) du3 



-, CO j 

1 X - r-* 



n J- 



L{C^{vi,V2,w)) dui- 



dt dv 



dv. 



Writing v = — sin 9 ■ vi — cos 9 ■ vj^ , we have 



du^' 



{-iy-^[sm9-^ + cos 



d 



d 



J-2 
1=0 



dvi dv^ 



COS 



Noting that v G L{C-^{vi,V2,w)) if and only if 6* G (0,7r — 9) we see 



j=2 •'■ 1=0 



I 



xHj.2-i{{vi,w))Hi{{vi,w))e-^\^^\"/^ 

X / sin^~2~'6'cos'0d6l. 
Jo 



Straightforward calculations show 

Kji{e) = {j-1) sm^^'^~^9cos^ede 
Jo 

ri-1 



/-B(j_i_,)/2,(/+i)/2(Vsin6'), if 6* > 7r/2, 

i-i)/2,{/+i)/; 
if (9<7r/2, 



p — I 

(-1)'^— (^(j-i-«)/2,(/+i)/2 - /-B(j_i_i)/2,(i+i)/2 (V sin 6* )) 



+ ^(j-l_i)/2,{«+l)/2 



where 



Jo 



is the incomplete beta function and B^^^^^^ = IBy-^^y^{l) is the beta function. 
Putting the above together, we have proved the following. 
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Lemma 5.4. The coefficient of / j\ in the power series expansion of 

lm.2{T{C{vi,V2,w),r)) 

is given by 

M]'''{C{VI,V2,W)) 



il-mvi,w))) 



'2tt 

{i-mvi,w))) 



2ir 



27r 



\/27r 



-y. 



1=0 



J -2 



^ ^ \ Kj^i{9)Hj_2-i{{vi,w))Hi{{vt,w))e 



-|hllV2 



j>2. 



It remains only to relate the above lemma to our original goal, that is, 
the EC densities of the field zi A Z2, which amounts to determining vi , V2 
and w for the cone K{u,p). We can take = (1,0) and V2 = {p, \/l — ) 
so that vi = (0, 1) and V2 = (a/1 — — p) and w = {u,u/^l + p). 
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